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3 I Abstract: 

We study the equation of state (EOS) of /3-stable dense matter and models of neutron 
stars in the relativistic mean field (RMF) theory with the isovector scalar mean field corre- 
sponding to the 5-meson [ao(980)]. A range of values of the 5-meson coupling compatible 
^ \ with the Bonn potentials is explored. Parameters of the model in the isovector sector are 
constrained to fit the nuclear symmetry energy, E s « 30MeV. We find that the quantity 
most sensitive to the 5-meson coupling is the proton fraction of neutron star matter. It 
increases significantly in the presence of the 5-field. The energy per baryon also increases 
but the effect is smaller. The EOS becomes slightly stiffer and the maximum neutron star 
Q-f mass increases for stronger 5-meson coupling. 

6 ' PACS: 21.65. -ff, 97.60.Jd 
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Nucleon matter in the RMF model with the 5-meson 



X 

The standard RMF model [1] of nuclear matter, frequently used in astrophysical cal- 
culations, involves mean fields of a, uj and p mesons. It does not include the contribution 
of the isovector scalar meson 8 [do (980)] although generally the density to which this field 
can couple does not vanish, < >^ 0. The contribution of the 5-meson field is not 

expected to be important for finite nuclei of small isospin-asymmetry, as the 5-meson mean 
field vanishes in symmetric nuclear matter. However, for strongly isospin-asymmetric mat- 
ter in neutron stars presence of the 5-field can influence the properties of dense matter. In 
Ref.[2] the RMF model was generalized to include the contribution of the 5-meson. Here 
we investigate consequences of such a generalized RMF theory for neutron stars. 
The dynamics of the RMF model is governed by the lagrangian 

L = L + Li nt , (1) 

where Lq = + L a + + L p + L§ is the free-field lagrangian and Li nt is the interaction 
term. The free-field lagrangians for nucleons and meson fields are: 
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H = i){id^ - m)ip, (2) 



where m is the bare nucleon mass, 




L a = ^(d^a-mla 2 ), 



(3) 



(4) 



(5) 



Here m CT , m Ul and m p are masses of respective mesons. Lagrangians (2)- (5) are the same 
as used in the standard RMF theory [1]. 

For the 5-field we use the simplest lagrangian of the massive isovector scalar field, 



where mj is the mass of the 5-meson. 

The coupling of the meson fields to nucleons is assumed to have the Yukawa form. For 
the u-field we use the cubic and quartic selfinteraction terms. The interaction lagrangian 
reads 



where U(a) is the potential energy term of the cx-field due to Boguta and Bodmer [3], 



In the RMF approximation, for a uniform nucleon matter only a few fermion densities 
are relevant. These include the baryon density, < V>7oV' > = n Bi the scalar density, 
< ijnp >, the isospin density, < ^qt^ >, and the scalar isospin density, < >. 
A selfconsistent description of the system is achieved by taking into account only those 
components of the meson fields which couple to the above densities with all remaining 
components vanishing. For normal nucleon matter the relevant components of mean meson 
fields are a, uq, p^\ and, for the 5-field, the isospin component 

It is a simple algebraic exercise to obtain the spectrum of nucleon energies in terms of 
the above components of meson fields and to construct all relevant nucleon densities. The 
single particle energies of protons and neutrons are 



L 5 = ±djd»6-±m 2 s 8 2 , 
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Lint = g*<rH> - 9ujU^^^ - -gppn^ril) + gsSiprip -U(a), 



(7) 




(8) 



Ep(N) (p) = 9^0 ± 2 9 P Po + y P 2 + m l 



■P(N) 



(9) 



where the proton and neutron effective mass is, respectively, 



mp = m - g a o - g^\ 



(10) 
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and 



m N = m - g a a + g 5 d^ s) . (11) 



The plus (minus) sign in the formula (9) refers to protons (neutrons). 
We obtain the field equations for mean meson fields in the form: 



mlu =g u n B , (12) 



where x = np/np is the proton fraction, 



ml& + — = g a (n s P + n s N ), (14) 



and 

In Eqs.(14) and (15) np and n s N is, respectively, proton and neutron scalar density, 



=g s (n s P -n s N ). (15) 



2 C ki m- 

< = 72^3 / d " k fur— -i ^ l = p ^^ ( 16 ) 

with proton and neutron effective mass, rap and ra^, given in Eqs.(lO) and (11). Mean 
meson fields as determined through Eqs.(12)-(15) depend on the baryon density np and 
the proton fraction x. 

The energy of the uniform matter consists of nucleon and meson contributions. The 
nucleon contribution is a sum of proton and neutron energies (9) up to their respective 
Fermi momenta. Mean meson field contributions are easily obtained from the lagrangians 
(4)-(6). With the mean meson fields given by Eqs.(12)-(15) the energy density of a uniform 
nucleon matter becomes 
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2 f^ p l f^ N I 1 

-^{J d 3 k^k 2 + m P + j d 3 k^k 2 + m 2 N ) + -Cln 2 B + 



+ ^[m - ^l^Lf + U{a) + ^d(2x - Ifnl + ^(m N - m P f . (17) 

In the spirit of the RMF theory [1] the parameters of the model are fit to repro- 
duce the empirical parameters of nuclear matter. In the isoscalar sector, the coupling 
parameters, C 2 = g 2 /m 2 ,C 2 = g^/m^, b = b/g 3 , and c = c/g*, are adjusted to fit the 
saturation properties of symmetric nuclear matter (x = 1/2), i.e. the saturation density 
no = 0.145/m~ 3 , the binding energy B = —16MeV per nucleon, and the compressibility 
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modulus K v R3 280MeV. The fourth parameter, e.g. c, can be used to control the stiff- 
ness of the equation of state of symmetric nuclear matter. Our choice of these coupling 
parameters is discussed below. 

The coupling parameters, C 2 = g 2 /m 2 , and Cf = g 2 /m 2 , in the isovector sector are 
constrained to fit the nuclear symmetry energy, E 8 = 31 ± AMeV [A]. This constraint 
gives C 2 as a function of Cf [2]. In the range of values of Cf considered below, which 
is compatible with the Bonn potentials [5], this function is well approximated by a linear 
relation [2] 

C 2 p =AC 2 s +B, (18) 

with positive coefficients, A > 0, and B > 0, depending weakly on the coupling parameters 
in the isoscalar sector. 

The requirement that the nuclear symmetry energy is reproduced in the presence of 
the 5-field is very important as the 5-field contribution is strongly attractive. As shown in 
Ref. [2] the 5-meson coupling provides a negative contribution to the symmetry energy that 
tends to cancel partially the p-meson contribution. Thus to preserve the nuclear symmetry 
energy a stronger p-meson coupling is needed. The formula (18) shows that the parameter 
C 2 has the lowest value C 2 = B for Cf = 0. For any coupling constant Cf > 0, the 
value of the p-coupling constant, C 2 , increases. Hence for pure neutron matter inclusion 
of the 5-meson results unavoidably in higher energy per particle at high densities, where 
contributions of vector mesons dominate. 

2. Equation of state and neutron stars 

As we mentioned above, there is one parameter, e.g. c, which can be used to label 
a family of EOS's in the RMF model. We wish to retain this freedom, as the true high 
density behaviour of the neutron star EOS is only weakly constrained at present [6]. In 
the following we shall mainly use two sets of coupling parameters which reproduce the 
saturation properties of nuclear matter but predict different stiffness of the EOS at higher 
densities. The EOS's corresponding to these two sets of parameters are referred to as soft 
and stiff. The soft EOS is specified by the parameters C 2 = 1.582/m 2 , C 2 = 1.019/m 2 , b = 
-0.7188, and c = 6.563. For the stiff EOS the parameters are C 2 = 11.25/m 2 , C 2 = 
6.483/m 2 , b = 0.003825, and c = 3.5 x 10" 6 . 

The above EOS's limit in a sense the family of EOS's in the RMF model. They are 
close to the softest and to the stiffest EOS in the RMF theory that are physically allowed. 
The stiff EOS corresponds to a very small value of the parameter c. Physical consistency 
of the RMF theory requires that c > 0. In Fig.l we show a plot of the parameter c as a 
function of C 2 . As one can see, c = for C 2 « 11.5/m 2 and thus any acceptable value 
of C 2 must be lower (Fig.l). In terms of the maximum neutron star mass, the EOS with 
c — >• is the stiffest acceptable EOS in the RMF model. We choose small but finite value 
c = 3.5 x 10 -6 with corresponding C 2 = 11.25/m 2 . 

The parameter c for the soft EOS has the highest allowed value, c = 6.563. The 
maximum neutron star mass corresponding to this EOS is about 1.44M . This is the 
mass of the heavier neutron star in the binary pulsar PSR B1913+16 [7], which has the 
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largest precisely measured value of a neutron star mass. Hence this EOS is the softest one 
still compatible with measured neutron star masses. 

We also show below some results for an intermediate EOS with the parameters C 2 = 
5.318/m 2 , C 2 = 2.31 /m 2 , b = -0.03952, and c = 0.4229. 

The coupling constant gs of the 5-nucleon interaction is a parameter in the one- 
boson-exchange fits of nucleon-nucleon scattering data. Its value is not, however, strongly 
constrained at present [5]. Our aim in this paper is to investigate the influence of this 
coupling parameter on the EOS of dense matter and on properties of neutron stars. To 
do so we adopt here a range of the 5-meson coupling compatible with the Bonn potentials 
[5], Cf < 4.4/m 2 . The maximum value of Cf we use exceeds the value corresponding to 
the Bonn potential C [5] which is Cf 2.5/m 2 . 

For a given 5-coupling, Cf , we obtain the p-meson coupling, C 2 , from Eq.(18) applying 
for the soft EOS the coefficients A = 0.63 and B = 5.0/m 2 . For the stiff EOS, the 
coefficients in Eq.(18) are A = 0.60 and B — 4.31/m 2 . 

To obtain the EOS of the neutron star matter we first calculate the proton fraction x 
of the charge-neutral /3-stable neutron star matter, which satisfies the condition 



f^N ~ VP = \J k% + m 2 N - yj k 2 p + m 2 P + ^C 2 p (l - 2x)n B = fi e = fi^, (19) 

where \i e = k e = (37r 2 n e ) 1 / 3 is the electron chemical potential and ^ = ^k 2 + m 2 is the 

muon chemical potential. Charge neutrality requires that n e + = np = xub- 

In Fig. 2 we show the proton fraction x as a function of density for a few values of 
Cf . One can notice that for both EOS's the proton fraction is substantially larger for 
indicated values of the 5-meson coupling, Cf, than for vanishing coupling, Cf = 0. It 
exceeds the critical value for the direct URCA process to dominate the cooling rate of 
neutron stars, which is xurca ~ 0.11, already at densities less than twice the nuclear 
saturation density. For the stiff EOS and for the 5-coupling corresponding to the Bonn 
potential C, Cf = 2.5/m 2 , the proton fraction is about 40% higher than for vanishing 5- 
meson coupling. For the soft EOS the proton fraction increases by about 20%, a somewhat 
weaker effect. 

The energy density of the /3-stable neutron star matter, e ns , is obtained as a sum of 
nucleon and lepton contributions, 

e ns = e nuc + ^kj + ^3 jf M d 3 k^k 2 + m 2 a7 (20) 

where nucleon contribution, e nuc , is given in Eq.(17). Next terms represent the energy 
density of the electron Fermi sea and the muon Fermi sea. The energy per baryon with the 
contribution of the 5-field for both soft and stiff EOS is shown in Fig. 3. For comparison, 
curves with no 5-field included are also shown. As one can see in Fig. 3 the energy per 
particle increases with 5-coupling Cf . The effect is stronger for the stiff EOS. For Cf = 
2.5/m 2 the energy per baryon is about 10% higher than for Cf = 0. For the soft EOS the 
energy per baryon increases by ~ 5%. 
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It is interesting to note that actually the nucleon contribution to the energy per baryon 
at a given ng decreases with increasing Cf for both EOS's. It is the lepton contribution 
which makes the total energy per particle higher. 

This behaviour can be easily understood. With increasing Cf the proton fraction 
of the neutron star matter in /3-equilibrium increases making the system less isospin- 
asymmetric. This in turn reduces the amplitude of the p-meson mean field, which 
provides a repulsive contribution. As a result, the nucleon energy per baryon is lowered in 
spite of the fact that the p-meson coupling parameter, C , is higher. However, increase of 
the proton fraction results in higher density of electrons and muons which thus contribute 
more to the total energy per particle. The energy per baryon of the /3-stable neutron 
star matter displays a similar behaviour as the one of pure neutron matter which becomes 
higher with increasing 5-meson coupling [2]. 

In Fig. 4 pressure as a function of mass density, P = P(e/c 2 ), a relation referred 
to as the EOS of neutron star matter, is shown. In this plot curves corresponding to 
Cf = AAfm 2 and to C$ = differ significantly for both soft and stiff EOS. Pressure at a 
given mass density is higher for Cf = AAfm 2 than for Cf = which indicates that the 
EOS becomes somewhat stiffer when the 5-meson contribution is present. 

The soft EOS and the stiff EOS calculated for Cf = 2.5/m 2 , a value corresponding 
to the Bonn potential C, are given in Table la and Table lb, respectively. As the proton 
fraction at high densities is rather high, especially for the stiff EOS, leptons (electrons and 
muons) make a sizable contribution to the total mass density. 

To gauge the influence of the 5-meson coupling on the EOS we have calculated models 
of neutron stars. The high density EOS calculated in the previous section was matched 
with the low density EOS due to Baym, Bethe, and Pethick [8] by constructing a proper 
phase transition. By making this construction we have found interface between crust 
matter, described by the EOS of Ref.[8], and liquid core matter, described by our RMF 
theory. Both EOS's in Table la and lb are given for liquid core matter starting from the 
core interface density. 

In Fig. 5 density profiles are shown for the canonical neutron star mass 1.4M©. For 
the soft EOS generally the radii are smaller than for the stiff EOS. The effect of the 8- 
meson contribution is more profound in case of the soft EOS. The radius increases from 
R « 10.6km, for Cf = 0, to R « 11.9km, for Cf = 2.5/m 2 , and to R « 12.5km, for 
Cf = AAfm 2 . For the stiff EOS the radius is R « 13.8km for Cf = 0, R « U.2km 
for Cf = 2.5/m 2 , and R xs 14.6km, for Cf = AAfm 2 . The central density decreases 
with increasing 5-coupling. In case of the soft EOS, the central density is n c « 8.5no, 
n c ~ 5.4no and n c ~ 4.4no for, respectively, Cf = 0, Cf = 2.5/m 2 and Cf = AAfm 2 . For 
the stiff EOS corresponding central densities are n c « 2.3n , n c 2.1n and n c 1.9n , 
respectively. 

It is interesting to study the core of 1.4M© neutron star models with the proton frac- 
tion exceeding the critical URCA value, xurca, where the direct URCA process dominates 
the cooling rate. With increasing Cf there are two effects influencing its size which tend to 
cancel one another: increase of the proton fraction at a given baryon density (Fig. 2) and 
decrease of the central density n c due to stiffening of the EOS. The former one tends to 
extend the core while the latter one tends to shrink it. As a result, the core size is rather 
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insensitive to changes of the 5-coupling Cf . It is determined primarily by the stiffness of 
the EOS for Cf = 0. Generally, the mass of the core where direct URCA process dominates 
is higher for the soft EOS than for the stiff EOS. 

The neutron star mass as a function of central density is displayed in Fig. 6 for the 
soft and the stiff EOS. We also show results for the intermediate EOS, with parameters 
given in Sect. 2. Maximum mass increases slightly with Cf. The 5-field plays a more 
important role for the soft EOS. The maximum neutron star mass is M max = 1.403Mq 
for Cf = 0. Since this value is less than the mass of the neutron star in the binary 
pulsar PSR B1913+16, which is 1.44M , this EOS is too soft to be realistic. However, 
with inclusion of the 5-field contribution the maximum neutron star mass increases to 
M max = 1.452M , for Cf = 2.5/m 2 , and to M max = 1.48M , for Cf = 4.4/m 2 . Inclusion 
of the 5-meson with Cf = 2.5/m 2 results in about 4% increase of M max that makes the 
soft EOS astrophysically acceptable. 

For the stiff EOS the maximum neutron star mass is M max = 2.275Mq, M max = 
2.3O9M and M max = 2.313M for, respectively, Cf = 0, Cf = 2.5/m 2 and Cf = 4.4/m 2 . 
In this case M max increases by ~ 1.5% for Cf = 2.5/m 2 . Central densities of the above 
maximum mass neutron stars are, respectively, n c ~ 6.2no, n c ~ 5.9no and n c ~ 5.7no 

3. Conclusions and discussion 

We have studied the influence of the 5-meson coupling on the EOS of neutron star 
matter in the RMF theory. When the isovector scalar field of the 5-meson is added to 
the standard RMF model, the nuclear symmetry energy has two contributions of opposite 
sign. The conventional p-meson contribution is positive, whereas the 5-meson contribution 
is negative [2]. This reflects the fact that in pure neutron matter the p-meson provides 
repulsion whereas the 5-meson produces additional attraction. It is thus obvious that 
in order to make physically relevant predictions of the EOS of neutron star matter the 
coupling parameters of the 6 and p mesons to nucleons should be constrained in such a 
way that the empirical value of the nuclear symmetry is preserved. 

We have found that the quantity most strongly affected by the presence of the 5-field 
is the proton fraction of the neutron star matter in /9-equilibrium. It increases rapidly with 
Cf . The effect is more profound for the stiff EOS, for which the proton fraction increases 
by ~ 40% for the value of Cf = 2.5/m 2 corresponding to the Bonn potential C. For the 
soft EOS the proton fraction increases in this case by ~ 20%. For higher values of the 
5-coupling parameter this increase is larger (Fig. 2). The URCA threshold concentration 
of protons, xjjrca ~ 0.11, occurs at lower densities. 

The 5-field contribution makes the energy per baryon of /3-stable neutron star matter 
higher by ~ 10% and ~ 5% for the stiff and soft EOS, respectively, and for the Bonn 
potential C value of Cf . Here the increase is due to larger leptonic contribution since the 
nucleon contribution actually decreases with Cf . The EOS with Cf ^ is slightly stiffer 
than for Cf = 0. Maximum mass of the neutron star also increases by a few percent. 

One should stress that the extension of the RMF model to include the isovector scalar 
meson 5 is the most natural one. In fact, the standard RMF model with no 5-field is not 
fully selfconsistent, as for isospin-asymmetric matter the density < tpr^ip >^ whereas 
the mean field with the same quantum numbers vanishes. The coupling constant g$ is not 
well determined at present [5] . However for values in the range compatible with the Bonn 
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potentials [5] the contribution of the 5-field in pure neutron matter at saturation density 
is quite strong [2]. If the p-coupling is unchanged, the neutron matter becomes selfbound 
for Cf ~ 1.0/m 2 . To avoid such an unphysical behaviour the p-coupling C 2 p should be 
readjusted to meet the requirement that the empirical value of the nuclear symmetry 
energy is reproduced. 
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Table la 

Soft equation of state, Cf = 2.5/m 2 



ns(cm 3 ) p(g/cm 3 ) P (dynes / cm 2 ) 
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Table lb 

Stiff equation of state, Cf = 2.5/m 2 
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Fig.3 

The energy per particle of the neutron star matter for the soft EOS and the stiff EOS. 
Curves labeled 1, 2, and 3 correspond, respectively, to Cf = 0, Cf = 2.5/m 2 , and C] = 
4.4/m 2 . 
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Pressure as a function of mass density for the soft EOS and the stiff EOS. Solid and dashed 
lines correspond, respectively, to Cf = and Cf = 4.4/m 2 . 
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Fig. 5 



The density profile of 1.4M neutron star for the soft EOS and for the stiff EOS. Solid, 
dashed and dotted-dashed lines correspond, respectively, to Cj = 0, Cf = 2.5/m 2 , and 
Cf = 4.4/m 2 . 
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The spectrum of neutron star masses for the soft EOS and the stiff EOS. Also, results 
for the intermediate EOS are shown. Solid, dashed and dotted-dashed lines correspond, 
respectively, to Cf = 0, Cf = 2.5/m 2 , and Cf = 4.4/m 2 . 
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